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Abstract 

Several examples of Jacobi matrices with an explicitly solvable spectral prob- 
lem are worked out in detail. In all discussed cases the spectrum is discrete 
and coincides with the set of zeros of a special function. Moreover, the compo- 
nents of corresponding eigenvectors are expressible in terms of special functions 
as well. Our approach is based on a recently developed formalism providing 
us with explicit expressions for the characteristic function and eigenvectors of 
Jacobi matrices. This is done under an assumption of a simple convergence con- 
dition on matrix entries. Among the treated special functions there are regular 
Coulomb wave functions, confluent hypergeometric functions, q-Bessel functions 
and q-confluent hypergeometric functions. In addition, in the case of q-Bessel 
functions, we derive several useful identities. 
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1 Introduction 



Special functions usually depend on a complex variable and an additional parameter 
called order. Typically, they obey a three-term recurrence relation with respect to the 
order. This is the basis of their relationship to Jacobi (tridiagonal) matrices. In more 
detail, the zeros of an appropriate special function are directly related to eigenvalues 
of a Jacobi matrix operator, and components of corresponding eigenvectors can be 
expressed in terms of special functions as well. One may also say that the characteristic 
function of the (infinite) matrix operator in question is written explicitly in terms of 
special functions. Particularly, Gard and Zakrajsek reported in [1] a matrix equation 
approach for numerical computation of the zeros of Bessel functions; on this point see 
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also [9J. In [8J, Ikebe then showed that the same approach was applicable, too, for 
determining the zeros of regular Coulomb wave functions. In practical computations, 
an infinite tridiagonal matrix should be truncated which raises a question of error 
estimates. Such an analysis has been carried out in [TO} IT2] . 

In [2], the authors initiated an approach to a class of Jacobi matrices with discrete 
spectra. The basic tool is a function $ depending on a countable number of variables. 
In more detail, we define $ '■ D — > C, 

oo oo oo oo 

#0) = 1 + 5^(-!) m ■■■ Yl x kl x kl+1 x k2 x k2+1 ...x km x km+1 , (1) 

m=l fci = l fc2=fci+2 fc m =fc m _i+2 

where the set D is formed by complex sequences x = {x k \ k K L l obeying 

oo 

Y \ X k%k+l\ < OO. (2) 
k=l 

For a finite number of variables we identify $(xi,X2, ■ ■ ■ , x n ) with $(x) where x = 
(xi,X2, • • • , x n , 0, 0, 0, ... ). By convention, we put $(0) = 1 where is the empty 
sequence. Notice that the domain D is not a linear space though £ 2 (N) C D. 

In the same paper, two examples are given of special functions expressed directly 
in terms of The first example is concerned with Bessel functions of the first kind. 
For w, v G C, v — N, one has 

Secondly, the formula 

00 jm(2m-l) 2m 

holds for t, w G C, \t\ < 1. Here o0i is the basic hypergeometric series (also called 
q-hypergeometric series) being defined by 

~ q k(k-i) 

and 

fe-i 

(a;q) k = ]J{l-aq j ), fc = 0,l,2,..., 

3=0 

is the g-Pochhammer symbol, see j5]. 

In [15], the approach is further developed and a construction in terms of $ of the 
characteristic function of certain Jacobi matrices is established. As an application, a 
series of examples of Jacobi matrices with explicitly expressible characteristic functions 
is described. The method works well for Jacobi matrices obeying a simple convergence 



2 



condition imposed on the matrix entries which is in principle dictated by condition 
(J2J) characterizing the domain of 

In the current paper, we present more interesting examples of Jacobi matrices 
whose spectrum coincides with the set of zeros of a particular special function. As a 
byproduct, we provide examples of sequences on which the function 5 can be evaluated 
explicitly. The paper is organized as follows. In Section [2] we recall from |141 ITS] 
some basic facts needed in the current paper. Section [3] is concerned with regular 
Coulomb wave functions. Here we reconsider the example due to Ikebe while using 
our formalism. In Section H] we deal with confluent hypergeometric functions. Here 
we go beyond the above mentioned convergence condition (see (j!4p below) which 
is violated in this example. Section is concerned with q-Bessel functions. This 
example is particular in that respect that the constructed second order difference 
operator is bilateral, i.e. it acts in £ 2 (Z) rather than in £ 2 (N). We first derive several 
useful properties of q-Bessel functions and then we use this knowledge to solve the 
spectral problem for the bilateral difference operator fully explicitly. Finally, another 
interplay between special functions, namely q-confluent hypergeometric functions, and 
an appropriate Jacobi matrix is demonstrated in Section [6j 



2 Preliminaries 

Let us recall from [T4"1IT5] some basic facts concerning the function $ and its properties 
and possible applications. First of all, quite crucial property of $ is the recurrence 
rule 

d({Xk}T=l) = 3tM£L 2 ) - X 1 X 2 #({^KL 3 ) ■ ( 5 ) 
In addition, $(xi, x 2 , ■ ■ ■ , Xk-i, Xk) = $(xk, Xk-i, . . . ,X2, X\). Furthermore, for x G D, 

lim $({x k }™ =n ) = 1 and lim $(xi,x 2 , ■ ■ ■ , x n ) = $(x). (6) 

Let us note that the definition of $ naturally extends to more general ranges of 
indices. For any sequence {x n }^= Nl , N±, ^ 6 ZU { — oo, +oo}, N\ < N2 + 1, (if Ni = 

N 2 + 1 G Z then the sequence is considered as empty) such that J2k=N\ \ x k%k+i\ < 00 
one defines 

00 

d[ixk}kl Nl ) = 1 + X!* -1 )" 1 Yl x kl x kl+1 x k2 x k2+1 . . . x km x km+1 

m=l kel(N 1 ,N 2 ,m) 

where 

I(N U N 2 , m) = {ke Z m ; kj + 2< k j+1 for 1 < j < m - 1, N t < k u k m < N 2 } . 
With this definition, one has the generalized recurrence rule 

$({ x k}kl Nl ) = ${{xk} n k=Nl ) $({ x k}kl n+1 ) - x n x n+1 ^({xkyiZ 1 ^) $({xk}kl n+2 ) 

(7) 
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provided n6Z satisfies N\ < n < N 2 . 

Let us denote by J an infinite Jacobi matrix of the form 



J 



f X\ W\ \ 

w 1 A 2 w 2 

W 2 A3 w 3 



(8) 



where {w n ; n G N} C C\{0} and {A n ; n G N} C C. In all examples treated in the 
current paper, the matrix J determines in a natural way a unique closed operator in 
£ 2 (N) (in other words, J min = J max ; see, for instance, [3J). If the matrix is real then 
the operator is self-adjoint. For the sake of simplicity of the notation the operator is 
again denoted by J. One notes, too, that all eigenvalues of J, if any, are simple since 
any solution {x k } of the formal eigenvalue equation 



A1X1 + wix 2 = zxi, Wk~ix k ^i + X k x k + w k x k +i = zx k for k > 2, 



(9) 



with z G C, is unambiguously determined by its first component x±. 

Let {7^} be any sequence fulfilling 7^7^+1 = w k , k G N. If J n is the principal nx n 
submatrix of J then 



det(J n — zl n 




In 



A, 



(10) 



The function $ can also be applied to bilateral difference equations. Suppose that 
sequences {w n }'^ = _ 00 and {( n }'^' = _ 00 are such that w n 7^ 0, ( n 7^ for all n, and 



E 



w 



CfcCfc+i 



< 00. 



Consider the difference equation 

w n x n+1 - ( n x n + = 0, n G Z. 



Define the sequence {V n } neZ by Vq = 1 and V n+ i = (w n /( n+ i)V n for all n. The 
sequence {7 n }nez is again defined by the rule 7„7 n +i = w n for all n G Z, and any 
choice of 71 7^ 0. Then the sequences {f n } n& and {g n } ne z> 



fn = V n $ 



21 00 



fc=n+l 



2 ^ n-1 



9r, 




(12) 



represent two solutions of the bilateral difference equation (rTT|) . With the usual defi- 
nition of the Wronskian, W(f, g) = w n (f n 9n+i — fn+ign), one has 



1 00 

' n=— 00 



(13) 
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For A = {An}^ let us denote Cg := C\ {A„; n G N}, and let der(A) stand for the 
set of all finite accumulation points of the sequence A. Further, for z G C \ der(A), let 
r(z) be the number of members of the sequence A coinciding with z (hence r(z) = 
for z G Cp). We assume everywhere that Cg 7^ 0. 

Suppose 



E 

n=l 



w 2 



(A n — 2o)(A n+ i — Zq) 



< 00 (14) 



for at least one z G Cq. Then (TT4"|) is true for all z G Cq [IS]. In particular, the 
following definitions make good sense. For k G Z + (Z + standing for nonnegative 
integers) and z G C \ der(A) put 



k 



&W : = .ta (u -, rW (n^M^ )■ (15) 



, =1 \iA,-u JZ=fc+1 
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Here one sets Wq ■— 1. Particularly, for 2 G Cq, one simply has 



Z=fc+1 



(this is in fact nothing but the solution /„ from f|12p restricted to nonnegative indices). 
All functions £fc(z), k G Z + , are holomorphic on Cq and extend to meromorphic 
functions on C \ der(A), with poles at the points z = X n , n G N, and with orders of 
the poles not exceeding r(z). This justifies definition (fTSj) . 

The sequence {£,k{z)} solves the second-order difference equation 

w k -ix k -i + (X k - z)x k + w k x k+ i = for k > 2. (17) 

In addition, (Ai — z)£i(z) + wi&iz) = provided £0(2) = 0. Proceeding this way one 
can show [T5| Section 3.3] that if £,o(z) does not vanish identically on Cq then 

spec(J) \ der(A) = {z G C \ der(A); £„(*) = 0}. (18) 

Moreover, if z G C \ der(A) is an eigenvalue of J then £(z) := (£i(z), £2(2), £3(2), ■ ■ •) 
is a corresponding eigenvector. If J is real and z G R fl Cg is an eigenvalue then 
HCWH 2 = £'o( z )£.i(, z )- Finally, let us remark that the Weyl m-function can be expressed 
as m(z) = fi(z)/fo(z)- 

Lemma 1. For p,r,£ G N, l<p<r + l<£, and any £-tuple of complex numbers 
%j, 1 < j < £> it holds true that 

XjX j+1 I d^X^Zl) ! S[{Xj} i j=r+ ^ 
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If p,r GN, 1 < p < r + 1, and a complex sequence {%j}°° =1 fulfills |||) then 



IOC 



\j=p-l 



Proof. Suppose {zj}°?__ 00 is any nonvanishing bilateral complex sequence. In [To^ 
Section 2] it is shown (under somewhat more general circumstances) that there exists 
an antisymmetric matrix Z(m, n), m,n G Z, such that 

3(m,n) = I ] [ — j $(z m +i, z m+2 , ■ ■ ■ , z n -i) 

\j=m+l 3 J 

for m < n, and 3(m, k)$(n, £) — 3(m, £)3(n, k) = 3(m, n)3(k, £) for all m, n,k,£ G Z. 
In particular, assuming that indices p, r, £ obey the restrictions from the lemma, 

3(0, r + l)3(p - 1, £ + 1) - 3(0, £ + l)3(p - 1, r + 1) = 3(0, p - l)3(r + 1, £ + 1). 

After obvious cancellations in this equation one can drop the assumption on nonvan- 
ishing sequences. The lemma readily follows. □ 

Lemma 2. Let x = {x n }^ =l be a nonvanishing complex sequence satisfying ([J|). Then 

F n :=$({xk}? =n ), n G N, (19) 

is the unique solution of the second order difference equation 

F n - F n+l + x n x n+l F n+2 = 0, nGN, (20) 

satisfying the boundary condition lim^oo F n — 1. 

Proof. The sequence defined in (lT9"j) fulfills all requirements, as stated in fl5j) 

and ([H])- It suffices to show that there exists another solution {G n } of f l2"U|) such that 
lim n ^ 00 G n = oo. If F\ = $(x) ^ then {G n } can be defined by G\ = and 

G « = f II —*—) ZU x *} n k=i) , forn > 2. (21) 

\fe=l X k X k+lJ 

If Fi = then necessarily F 2 ^ since otherwise (120]) would imply F n = for all rz 
which is impossible. Hence in that case one can shift the index by 1, i.e. one can put 

G 2 = 0, 

'n-2 



G »= ( II ) ^(WJS) » forn > 3, 

(and G*! = — Xix 2 ). In any case, F n is the minimal solution of (T2"0"j) . see [5J. □ 
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Remark 3. If $(x) = then $(xi, x%, ■ ■ ■ , x n ) tends to as n — > oo quite rapidly, more 
precisely, 

$(x 1 ,x 2 ,...,x n+1 ) = o^JJ^x k x k +ij, as n -> oo. (22) 

In fact, if 5(2?) = then F2 7^ and the solutions {F n } and {G n } defined in (fT9]) and 
(|2Tj) . respectively, are linearly dependent, F n = F 2 G n) Vn. Sending n to infinity one 

gets 

l = F 2 lim fn^— ]?(W1U)= Mm ( ]1 W({x fe }Li) ^(WS) ■ 

Now, Lemma [1] provides us with the identity 

n 

m*k}u)9({*k}Ei)-d({xk}Ei)m*k}u = i[x k xk + x, 

k =i 

and so one arrives at the equation 



fc = 2 ) = 0. 



Since #({x fc }£° =2 ) ^ this shows (F22 



3 Coulomb wave functions 

For x > 1, y G R, put 



A(x,y) 



?7 



and 



(x — l)x ' 



w(x,y) 



1 

X 



x 2 + y s 



y/2x~=ir(l(x + 1)) 



4a; 2 - 1 ' 
r(|(x + ^ + l)) 



r(|(x + ^)) 



Then 7(2?, y)^(x + l,y) = w(x,y). For > 0, v G R, consider the Jacobi matrix 
J = J(/i, 1/) of the form ([8]), with 



A fc = A(/i + k, v), w k = w(p + k,v), k — 1, 2, 3, 



(23) 



Similarly, 7^ = 7(/i + /c, z/). Clearly, the matrix J(/x, z/) represents a Hermitian Hilbert- 
Schmidt operator in £ 2 (N). Moreover, the convergence condition (j!4jl is satisfied for 
any z £ C\{0} such that z 7^ X k , Vfc G N. 

Recall the definition of regular Coulomb wave functions [H Eq. 14.1.3] 



F L ( V ,p) = 2 L e~^ 2 



\F(L + l + ir,)\ ^ 



T(2L + 2) 



p L+L e- lp \Fi(L + l-i7 1 ;2L + 2; 2ip), (24) 



valid for L G Z + , r\ G IR, p > 0. Let us remark that, though not obvious from its 
form, the values of the regular Coulomb wave function in the indicated range are real. 
But nothing prevents us to extend, by analyticity, the Coulomb wave function to the 
values L > — 1 and p G C (assuming that a proper branch of p L+l has been chosen). 

As observed in [8J, the eigenvalue equation for J(fi, v) may be written in the 
form F^_i(— i/, = 0. Moreover, if z ^ is an eigenvalue of J(/x, v) then the 
components v n (z), n G N, of a corresponding eigenvector t> (z) are proportional to 
y/2fi + 2n — 1 Ffj, +n -i(— u, z^ 1 ). Thus, using definition (1241) . one can write 

spec( J(ji, u))\{0} = {C 1 ; e"* 1 F 1 (ji + iu; 2/x; 2iQ = 0} (25) 

and 

v n {C X ) = v /2 / x + 2n-l |r ^ + ^+^ )l (2C)- 1 ^(ai + n + ii/; 2// + 2n; 2iQ. 

^ U (26) 
Here we wish to shortly reconsider this example while using our formalism. 

Proposition 4. Under the above assumptions (see 

I^-tV ) (27) 



h - C" 1 



k=l 



r^roz + i) e i^(^+^ 2 ^ 2 <)- 

Proof. Observe that the convergence condition (jl4p is satisfied in this example. For 
n G N put 



/l,n = $ 



Ik 2 



k=n 



and let f2, n be equal to the RHS of fl27|) where we replace /i by + n — 1. According 
to (J^J), the sequence {/i, n } obeys the recurrence rule 



where 

X(x) 



fi,n - + X{n + n)fi, n+2 = 0, n G N, (28) 



(A(x,z/)-C- 1 )(A(x + l,z/)-C- 1 ) 
(a; 2 -l)(x 2 + z/ 2 )C 2 



(4x 2 - 1) {{x - l)x - i/C) (x(x + 1) - i/C) 
Next one can apply the identity 



for x > 1. 



rn/ n , on 6 2 - 26 + (2a - 6);z . , . a(b-a)z 2 _. 1 , n , 
iF^o-l; 6-2; — ^ 1 F 1 (a; 6; - t^j^ i*i(a + l; 6 + 2; = 0, 

(29) 
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as it follows from jU §13.4], to verify that the sequence {f'2,n} obeys (I28|) as well. 
Notice that, if rewritten in terms of Coulomb wave functions, (129|) amounts to the 
recurrence rule pQ Eq. 14.2.3] 

Ly/(L + iy + ^u L+1 - (2L + l)L+ L(L ^ +1) ^ L + (L + 1) V^T^ML-l = 0, 

where = Fl(t],p). 

To evaluate the limit of / 2i „, as n -> oo, one may notice that 

lim iFx(a + n;b + «n; 2) = e z/K 

for k^O, and apply the Stirling formula. Alternatively, avoiding the Stirling formula, 
the limit is also obvious from the identity [7, Eq. 8.325(1)] 

n (1 1 2 \ rwrfa+D m 

li V (v + k)( v + k + i)J rd + j-ivr^rg+j+ivni)' v ' 

In any case, lim^oo / 2 ,n = 1 and so, in virtue of Lemma [2J / l n = / 2 n , Vrt. In 
particular, for n = 1 one gets f[2"Tj) . □ 

Proof of formulas < fl?5]) and (R^j . As recalled in Section [2] (see f llSp ). 2 = £ _1 7^ is 
an eigenvalue of J(fi,i>) if and only if £,o{ z ) = which means nothing but (I25p . In 
that case the components £ n (z), n G N, of a corresponding eigenvector can be chosen 
as described in (I15D. Note that 



n-1 



T(yU + n + zV) 



r(/x + i + ^) 



r(2 At + i) 



T(2/i + 2n) 



J^[w fc = 2 n -V(2^ + l)(2/^ + 2n- 1) 

fc=i 

and that (130]) means in fact the equality 

yr 1 r(| + /i - \ yi + 4z/z) r(| + /i + |Vj + 4^) 

1 1 1 - + fc, i/) 2r ~ r(/x)r(/x + 1) ' 

Using these equations and omitting a constant factor one finally arrives at formula 
flUD. □ 

4 Confluent hypergeometric functions 

First, let us show an identity. 
Proposition 5. The equation 



fc=i/ 



r(x + 7 + n) /] v / 2^r(|( 7 -a + fc + 1)) 



r(x + 7 ) ^[(x + 7 + A;- l)r(i( 7 -a + A;)) 

p(^y + n) 

— — - — iFi(a; 7; 2) iFi(a - 7 - ra; 1 - 7 - n; x) (31) 
r (7) 

Yt/ — n/ -I- n -I- 1 

Fi(a - 7 + 1; 2 - 7; 2) iFi(a; 7 + 71 + l;x), 



r( 7 -l)r( 7 -a + n+l) „ +1 



r(7-a)r( 7 + r2 + l) 



is valid for a,j,x G C and n G Z + (if considering the both sides as meromorphic 
functions). 



Remark 6. For instance, as a particular case of f l3T|) one gets, for n = 0, 

fry Ql\X 

iFi(a;7;x) iF x {a -7; 1 -7; x) — iFi(a; 7 + 1; x) iFi(a~7 + l; 2-7; x) = e x . 

7(7 - 1) 

(32) 

Proof. For a, 7 and x fixed and n£Z, put 

<^> n = tt7 ~~ — r iFi(a;n + 7;x), ^ n = — — — r U{a, n + 7, x). 

r(n + 7) r(n + 7 — a) 

Then {tp n } and {ip n } obey the second-order difference equation [1, Eqs. 13.4.2,13.4.16] 

(n + 7 — a)xu n+ i — {n + 7 + x — l)u n + w„,_i = 0, n G Z. (33) 

Note also that 

r(7 + i) 

(a - 7) iFi(a;7+ l;x)?7(a,7,x) + 7i-Fi(a;7;x)[/(a,7 + l,x) = — — — — x 1 e x 

T(a) 

(as it follows, for example, from equations 13.4.12 and 13.4.25 combined with 13.1.22 
in pQ). Whence 

r(a)r(l - a + 7) 

and so the solutions tp n , ip n are linearly independent except of the cases —a G Z + and 
a - 7 G N. 

The difference equation (133]) can be symmetrized using the substitution 

x -n 

T(7 — a + n) 

Then w n v n+1 - ( n v n + w n _iw„_ 1 = where 

x~ n (x + 7 + n — l)x~" 

n r(7 — a + n) ' n r(7 — a + n) 

For a solution of the equation 7„7 n +i = w n , Vn, one can take 



1 nil 



7„ = 2±x 2 + 4 . 



r(i( 7 -a + n+ 1)) 



T(7 - a + n)r(|(7 — a + ri)) 
Referring to another solution, namely 

-1 / f 2 ~\ n—l 

1 */ \lk 



9n = f4 h with n G N, 
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using otherwise the same notation as in (112)) . one concludes that there exist constants 
A and B such that 



r(7 - a + n + 1) \ I (x + 7 + k - 1) r(§(7 -a + k)) 



for all n G Z + . A and i? can be determined from the values for n = —1,0 (putting 
$({xk}k-Li) = 0) as dictated by the recurrence rule (j7j) provided the admissible values 
are extended to N± = N2 = 1, w = 0). After some manipulations one gets 

IY7 — ct) 

iFi(a;7;a;)C/(a,7 + n + 



r(7 — a + n + 1) 

r(7) 



r( 7 + n + l) 



f/(a,7,x) 7 + n + 1; x) 



r( 7 )r( 7 -«)r(x + 7 + n) ^_ 7 _ n /j v^r(§( 7 -« + * + !)) 



fc=i/ 



r(a)r(7 - a + n + l)r(z + 7) u \J (x + 7 + k - l)r(±( 7 - a + Jfe)) 
Recall that 

U(a, b, x) = J {1 ~ b \, 1F1 (a; 6; x) + ^ ^(a - 6 + 1; 2 - 6; x), (34) 

T(a — o + l) r(a) 

whence (1311. □ 



Remark 7. Let us point out two particular cases of ( 13T|) . Putting a = one gets the 
identity 



v^r(i( 7 + fc + i)) \ r(x + 7 ) A r( 7 + n-j) , 
* \j (x + 7 + fc - i)r(|( 7 + k)) j k J r(x + 7 + n) r( 7 ) ' 1 ' 

and for a = — 1 one obtains 

v / 2^r(|( 7 + fc + 2)) 



(x + 7 + A;-l)r(i(7 + fc + l)) 
T(x + 7) v^T(7 + n-j) 



E rrIV (7-j(n-j))^- (36) 



r(x + 7 + n) ^ r( 7 + i) 

Let us sketch a derivation of f|5B|) . equation (|35|) is simpler. Substitute —1 for a 
and 7 + n for 7 in (13"2"|) , and put 

= r (7 + n ~ ^ a.-n lFl (- 7 - n; 2 - 7; x). 
7 — x + n 



R R r( 7 + n + l) n-i, 
-Dn+i — £>n — ~ — — x e 



Then 

R_ , , - R._ = 

(7 — a; + n) (7 — a; + n + 1) 
11 



Whence 

%i - -do + e > t — tt-7 —t- x 

(y-x + j - l)(y-x + j) 

which means nothing but 

T(7 + n)(7 — x) 7 — n — 1; 1 — 7 — n; x) 

- T( T - 1)( 7 - x + n + 1) x™ +1 ^(-7; 2 - 7; x) (37) 



( 7 -x + n + l)( 7 -x)e a; V- 

( 7 -x + j-1)(7-x + j) 



Set a = — 1 in fl3T|) and notice that i-Pi( — 1; b] x) = 1—x/b. After some simplifications, 
a combination of thus obtained identity with (j37j) gives fl36|) . 

As an application of (l3Tj) consider the Jacobi matrix operator J(a, /?, 7) depending 
on parameters a, (3, 7, with (3 > 0, 7 > and a + (3 > 0, as introduced in (jHJ) where 

we put 

A fc = 7/c, w fc = \fa + fc = 1, 2, 3, — (38) 

For the sequence 7^ (fulfilling 7^7^+1 = Wk) one can take 

V = ,Wi(H + fc + iY)/rfi(£ + * 



2 V./3 /// V 2 V/5 

Regarding the diagonal of J(a,/3, 7) as an unperturbed part and the off-diagonal 
elements as a perturbation one immediately realizes that the matrix J(<x/3,7) deter- 
mines a unique semibounded self-adjoint operator in £ 2 (N). Moreover, the Weyl theo- 
rem about invariance of the essential spectrum tells us that its spectrum is discrete and 
simple. Our goal here is to show that one can explicitly construct a "characteristic" 
function of this operator in terms of confluent hypergeometric functions. 

Proposition 8. The spectrum of J(a, /3, 7) defined in and ( Tffgp coincides with the 
set of zeros of the function 



^ a , A7i , )=l , v |i-|-£-£ i i-£-£ i £)/r(i-A-£). (39) 



Moreover, if z is an eigenvalue then the components of a corresponding eigenvector v 
can be chosen as 



Vk 



r(g + * 



1/2 



D'^V , K \ ' M f.(i-g-4- £ ;i-4- £ + *;4 



r^i-A_£ + j^j V. /? 7 2 7 7 2 7 7 2 , 
fceN. (40) 
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Remark 9. (i) In principle it would be sufficient to consider the case 7 = 1; observe 
that 

Fj(a, B, 7; z) = Fj ( — , 1; - 

V7 2 7 2 1, 

Thus for 7 = 1 we get a simpler expression, 

Fj(a, 8, 1; z) = x Fi ( 1 - ^ - - *; 1 - /3 - z; s\ j V{1 - 8 - z). 

(ii) Notice that the convergence condition (Il4p is violated in this example. 

Before the proof we consider analogous results for finite matrices. Let J n (a,/3,7) 
be the principal n x n submatrix of J[a,8,^f). The characteristic polynomial Fj n (z) 
of J n (a,8,j) can be expressed in terms of confluent hypergeometric functions, too. 
According to ( JTUj) . 

2g r(|(g + fc + i 

" (*-*MK? + *))/_,< 

Applying (1311) one arrives at the expression 

7;*) 

rfn + i-4 





7 ™ e -^| r ^ ljFl A - ^ - -g- - £; 1 - -g- - £; A 

' JL _ £^ V /3 7 2 7' 7 2 7' 7 2 



r 1 - 4 - £ 

7 i 7 



_, / a (3 z 8 

x jFi -n - -; -n + — + -; — 

\ p 7 Z 7 7^ 

„ +i r („ + i + g )r(-^- f ) , «.,,£,*.£ 
W r ( ? )r( B + 2-^-i) 1 7 + y + ?7 2 



" _ A _ £■ _L9_il_£-A 

p 7 7 7 7 7 



Eigenvectors can be explicitly expressed as well. If z is an eigenvalue of J n (a, 8, 7) 
then formula (TT6|) admits adaptation to this situation giving the expression for the 
components of a corresponding eigenvector, 



Cr = (- 1 ) 111^)1 II (Ai-^^Mr^- |, * = l,2,...,n. 

A i - z ) j= k+i 

Notice that £^ makes sense also for k — n + 1, and in that case its value is 0. Using 



31]) and omitting a redundant constant factor one arrives after some straightforward 
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computation at the formula for an eigenvector of J n (a, 7): 

(n) = (_l)kpk/2^-k 1 i \P J \ 7 7 J 



V k 



r[i-£-* + * 



a /9 2 ., /9 z , , A P / a /? . z 
X 1F1 ( 1 - ; 1 + k; — 1F1 -- - n; — + - - n; — 

p 7^ 7 7 7 7 / VP 7 7 7 







7 ^ r(2-£-* + /, 

XiF 1 ( r i-^-A;;l + 4 + --*;4 

p 7 Z 7 7^ 

a /3 z z \ \ 
x 1 F 1 l---±L--;2--%-- + n;4 , 1< A; < n. 



Remark 10. Formula (HO]) can be derived informally using a limit procedure. Suppose 
z is an eigenvalue of the infinite Jacobi matrix J (a, 0, 7). For fceN fixed, considering 
the asymptotic behavior of v% as n — > 00 one expects that the leading term may 
give the component v k of an eigenvector corresponding to the eigenvalue z. Omitting 
some constant factors one actually arrives in this way at ()40|) . But having in hand the 
explicit expressions (l39j) and (j40j) it is straightforward to verify directly that the former 
one represents a characteristic function while the latter one describes an eigenvector. 

Proof of Proposition Observe first that for k = the RHS of ( HO]) equals, up to a 
constant factor, to the announced characteristic function ( I39p . If z solves the equation 
Vq = then one can make use of the identity [JJ Eq. 13.4.2] 

b(b - 1) x Fi(a; 6 - l;x) + 6(1 - b - x) iF^a; 6; x) + (6 - a)x 1 F l (a; 6 + 1; 2) = 

to verify that t> G £ 2 (N) actually fulfills the eigenvalue equation (J9j. Note that the 
Stirling formula tells us that 

:-i)* , . . 4! /j e \"7, . <\ 



On the other hand, whatever the complex number z is, the sequence v k , k G N, solves 
the second-order difference equation (]T7]) . and in that case it is even true that 

U>OVo + (-^1 — z)Vi + WiV 2 = 0. 

Let g k , fc G N, be any other independent solution of (fT7|) . Since the Wronskian 

w k (v k g k+1 - v k+1 g k ) = const ^ 

does not depend on k, and clearly lim^oo w k v k = lim^oo w k v k+ i = 0, the sequence 
g k cannot be bounded in any neighborhood of infinity. Hence, up to a multiplier, {v k } 
is the only square summable solution of (ITT]) . One concludes that z is an eigenvalue 
of J (a, (3, 7) if and only if wqVq = (which covers also the case a = 0). □ 
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Remark 11. A second independent solution of ()17|) can be found explicitly. For exam- 
ple, this is the sequence 

gt = ( _ 1)r / V 'r(^ + kY 1,2 u( i _ 2 _ 4 _ i,i _ £ _ i + *, 4), fc 6 N , 

\P J \ P 7 2 7 7 7 7 / 

as it follows from the identity [TJ Eq. 13.4.16] 

(6 - a - l)C/(a, 6 - 1, x) + (1 - 6 - x)U(a, b, x) + xU(a, b + l,x) = 0. 
But using once more relation (|34jl one may find as a more convenient the solution 



9k 



P , z k 

/-y2 



r(f + A;)r(l + 4 + £-A; 



Remark 12. Let us point out that for a = one gets a nontrivial example of an 
unbounded Jacobi matrix operator whose spectrum is known fully explicitly. In that 
case 

A fc = jk, w k = \ffik, k = 1,2,3, ... , 

and 

P z 



Fj(0,(3, r ,z) = e^ 2 / /r^l 



7 2 7 



Hence 

spec J(0, P, 7 ) = | + 7 j; J = 1, 2, 3, . . .| . 

Remark 13. Finally we remark that another particular case of interest is achieved in 
the formal limit ft — » 0. Set a = w 2 for some w > 0. Since [U Eq. 13.3.2] 

lim ^Ja-b---) =z^ 2 T{b)J b ^{2V~z) 

one finds that 

fw\ zh f2w\ 
lim Fj( W 2 , /3, T ,z) = ( - J J- Zh ( — ). 



/3^o \7/ V 7 / 

It is known for quite a long time jU [H] that actually 



spec J(w 2 , 0, 7) = < z G C; J. 



-2/7 



2w 

7 
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5 Q-Bessel functions 

5.1 Some properties of q-Bessel functions 

Here we aim to explore a q-analogue to the following well known property of Bessel 
functions. Consider the eigenvalue problem 

wxk-i — kxk + wxk+i = vxk, k G Z, 

for a second order difference operator acting in £ 2 (Z) and depending on a parameter 
w > 0. If v ^ Z then one can take {J u+k (2w)} and {{—\) k J _ u _ k {2w)} for two 
independent solutions of the formal eigenvalue equation while for v G Z this may be 
the couple {J u+ k{2w)} and {Y u+ k(2w)} . Taking into account the asymptotic behavior 
of Bessel functions for large orders (see [U Eqs. 9.3.1, 9.3.2]) one finds that a square 
summable solution exists if and only if v G Z. Then x^ = J u+ k{2w), k G Z, is such a 
solution and is unique up to a constant multiplier. Since 

oo 

E J *w 2 = X > ( 41 ) 

fc=— oo 

thus obtained eigenbasis u v = {?V,fc}&l_oo) ^ £ Z, with tv^ = J v+ k(2w), is even 
orthonormal. One observes that the spectrum of the difference operator is stable and 
equals Z independently of the parameter w. 

Hereafter we assume < q < 1. Recall the second definition of the q-Bessel 
function introduced by Jackson [TT] (for some basic information and references one 
can also consult [5j), 



4 2 \x;q) 



(q u+1 ;q)oo rxy , ( „ +1 <f 



(g;g)oo V2, 

Here we prefer a slight modification of the second q-Bessel function, obtained just by 
some rescaling, and define 



_„v+3/2 ±_ 



Ux; q) := ^AjfV /4 *l q) = q^^ ( ^j g) °° (|)" g' +1 ; g, -g 4 

? ' 9 °° (42) 
With our definition we have the following property. 

Lemma 14. For every nGN, 

U(x;q) = (-irUx;q). (43) 

Proof. One can readily verify that 

„n 2 / 2 7/, 2n 

urn ii — q ) <M, g , g, — g w j — — - — r — -; — — <M, g , g, — g w 

(g; g)n-i(g; g)n 

and 

The lemma is an immediate consequence. □ 
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Proposition 15. For < q < 1, w, v G C, q v £ q z+ , one has 



y (| ^pqw )J = °M;^-f +1/2 ™ 2 )- (44) 

Remark 16. If rewritten in terms of q-Bessel functions, (I44p becomes a q-analogue of 
©. Explicitly, 



where [5] 

\ x \ = C'^l rjx) = (i - q y~*. 



Lemma 17. For v G C, q v q z+ , and all s6N, 

(<T+ fc ; q) s+1 ~ (1 - q>) (</": </) 



fc=0 

Proof. One can proceed by mathematical induction in s. The identity 

gsfc q(s~l)k ( \ 1 



(^+ fc ;g) s+ i g"(l-g s ) V(g" +fe ;g) s (g ! ^ +fe+1 ;g) i 
can be used to verify both the case s — 1 and the induction step s — >■ s + 1. □ 

Proof of Proposition [731 One possibility how to prove (|44j) is based on Lemma [2J The 
proof presented below relies, however, on explicit evaluation of the involved sums. For 
v G C, q v i q z , k G Z, put 



Pk 



q (u+k)/2 



1 - <f + fc ' 

Then ([45]) immediately implies that, for nfZ and s G N, 

00 s(i/+n+l)/2 

E q {s-l){p+k)/2 _ _ _ _ _^ p n p n+ i . . . p n+s -i. 

1 — q s 

k=n ^ 

This equation in turn can be used in the induction step on m to show that, for m G N, 
n E Z, 

oo oo oo 

y~] ■ ■ ■ 5j PkiPkx+iPk 2 Pk 2 +i ■ ■ ■ Pk m pk m +i 

k 1 =nk 2 =k 1 +2 fc m =fc m _i+2 



^m(3m+l)/4 ~m(i/+n-l)/2 



In particular, for n = 1 one gets 



' PnPn+l ■ ■ ■ Pn+m—1- 



00 00 00 ^m(2m+l)/2+™ 

E E '■■ S p kl p kl+1 p k2 p k2+1 • . . p fcm p fcm +i = 7— r 7 — I , , m G N. 

fe 1= lfc 2 =fel+2 fc m =fc m _!+2 W '^ mW ' yM 

Now, in order to evaluate ^{wp^}^), it suffices to apply the very definition ([1]). □ 
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The q-hypergeometric function is readily seen to satisfy the recurrence rule 

o0i(; <t\ q, z) - o 0i(; q u+1 ; q, qz) - jz -r-. — tt <M; q v+2 \ q, q 2 z) = o. 

(1 — q v ) (1 — g^ 1 ) 

Consequently, 

w) u {2w; q) - (q-^/ 2 - g^+D/ 2 ) j u+1 (2w; q) + w) v+2 {2w; q) = 0. 

This is in agreement with (TT2I) if applied to the bilateral second order difference equa- 
tion 

wVj - (q^+ n )/ 2 - g ("+*)/2) Xn + =0, fiGZ. (46) 

Suppose q u ^ q z . Then the two solutions described in ()12j) in this case give 

fn = q- v ^W ( { ™\ w-") v ^(2w;q), (47) 

{q u+ ,q)oc 

9n = (-ir+V^ +1)/4 7%^y ^j-,- n (2 W ;g), n G Z. (48) 

Let us show that they are generically independent. For the proof we need the identity 
13 §1,3] 

^ q k[k-l)/2 



fc=0 



Lemma 18. For v G C, q u ^ g z ; the Wronskian of the solutions of (Jty ), {f n } and 
{g„} defined in fcffy and (fify , respectively, fulfills 

Proof. The first equality in (1491) is nothing but (TT3j) . Further, in virtue of (T4*4|) . the 
second member in (149)) equals 

oo k(k— 1) 

M ^°° ~ (q;q)k(q- M ;q)k K ' 

00 „fc 2 /2 

Ey 2fc / 1/2 2 \ 
te „fe^ = <-*'«">-■ 

The lemma follows. □ 
At the same time, W(/, 5) equals 

-i/(H-l)/2 / \ 2 

L- +1 -g)oo(g - g)I (U2u;; g) ^- l{2w] q) + Ul(2w; q) ^ {2w; q)) ■ 

This implies the following result. 
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Proposition 19. For w G C one has 

j v (2w; q) j^_i(2w; q) + j v+ i(2w; q) )- v (2w; q) 



(g;g)oo 2 w; 

and, rewriting / T577j) m terms of q-hypergeometric functions, 



(50) 



o<M;g" +1 ;g, -g^ + M o<M; g "«) 
— it o0i(; g^ 2 ; g, -g^) o<M; g^ +1 ; g, -g" 1 ^) (51) 



(1 - q u ) (1 - q u+ 

Remark 20. Let us examine the limit q — > 1— applied to flSTJj) while replacing w by 
(1 — q)w. One knows that [5] 

lim )„((1 - q)z- q) = J v (z), lim (1 - qf- = !». 

Thus one finds that the limiting equation coincides with the well known identity 

J v (2w) J- v -!(2w) + J v+l (2w)J_ v (2w) ' MU( T?/) 



iur(i/ + i)r(— i/) 7rif 

It is desirable to have some basic information about the asymptotic behavior of 
q-Bessel functions for large orders. It is straightforward to see that 

- )v ( X ;q) = q ^+W_L_ (*Y(i + (f)) as Rez/^+oo. (52) 
{q;q)oo ^2/ 

The asymptotic behavior at — oo is described as follows. 
Lemma 21. For a, w G C, q a ^ q z , one has 

lim siYi{-Kv)q v{u+1)/A w- u ) u {2w-q) 

vE— a— N 

= - sin(™) g^-V* ^ g^(g^ % tf^l *)~ . (53 ) 

(g; g)oo 

Proof. Put i/ = — er — n where nGN. Using flH]) and (J7|) one can write 

o<M; <?-"-";<?, - g - CT - n+1 / 2 w 2 ) 

= ^({ g( CT + fc )/2 _ q -{°+k)/2 } J ^({ g(^-fe)/2 _ q-(°-l°)/2 J 

+ 



(g<7/2 _ g -a/2^ g (l~a)/2 _ g -(l-a)/2^ 

x s /„ xl „, ._c^_i_m/o r I 5 



q(v+k)/2 _ q -{a+k)/2 j J q(<r-k)/2 _ g -(a-k)/2 } 
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Applying the limit n — > oo one obtains 

lim oM;q~ a - n ;q,-q- a ~ n+1/2 ™ 2 ) 

n— »oo 

= oU; <f\ 9, -q^'w 2 ) o 0i(; g 1 ^; q, -q^ +3/2 w 2 ) 
+ 



w; 2 



_ 9 -«r/2)(g(l-(0/2 - ? -(l-«r)/2) 

x g 1+ - g, -g-+ 3 / 2 u, 2 ) O X (; g 2 ~- g, -g^V) 

To get the last equality we have used (1ST]) . Notice also that 

Jim n+l)/2 ( g — g)oo = ^«r-l)/2 g) 00 (g 1 "- g)^. 

n— >oo 

The limit (j55|) then readily follows. □ 

Finally we establish an identity which can be viewed as a q-analogue to flUJ . 
Proposition 22. For < q < 1 and wfC one /ias 

OO oo 

£ g- fc / 2 j fc (2u;; g) 2 = j (2u;; g) 2 + £ (g fc / 2 + g^ 2 ) j fc (2«,; g) 2 = (-g 1 ^ 2 ; g )oo . 

k=—oo k=l 

(54) 

Equivalently, if rewritten in terms of q-Bessel functions, 

oo 

4 2 \2w, q) 2 + J2 (q k/2 + q- k/2 ) q k2/2 4 2 \^ q) 2 = (-^ 2 ; <?)«, 
fc=i 

Proo/. In [JJJ (1.20)] it is shown that 

J '' ,{2lr:,l)2 - A^ +1 ^)-V^ 3 2 (g^,-g^,-g^;g^,g 2 ^;g,- W ; 2 ) ! 



(-u> a ;g)oo V (?;g)oo 

and this can be rewritten as 

O0i(; q u+1 ; q, ~q u+1 x) 2 = (-*; g)^ Md+l-if+l, -q u+1 ; q u+ \ q 2u+1 ; q, -x). 
Hence (113]) is equivalent to 

302(g 1/2 , -g 1/2 , -q; q, g; g, -x) 

00 ^- fc / 2 4- ^/ 2> ) « fe 7 2 

+ E~ — 302(g A:+ i-? fc+ i-g fc+1 ;g fc+1 ,g 2fc+1 ;g ) -a;)x fe = l. 
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Looking at the power expansion in x one gets, equivalently, a countable system of 
equations, for n — 1, 2, 3, . . ., 

(q 1/2 ;q)n(-q 1/2 ;q)n(-q]q)n 
(q; q)n 

, 1 )k («T* /2 + d k ' 2 ) a* 2 ' 2 (^ +1/2 ; q)n- k (-q k+1/2 ; q)nM-q k+l ; q)n-k = Q 

fcl ' (q;q)k 2 (q;q)n-k(q k+1 ;q)n- k (q 2k+1 ;q)n^ 

The equations can be brought to the form 

i 1)k q k(k - 1)/2 (i + q k ) _ 



(<?;<?)n 2 ^ (g;g)n+fc(g;g)n-fc 



or, more conveniently, 



2" -j(2n-j+l)/2 

y^(_i)i_« = o. 

^ (q;q)2n-j(q;q)j 



This is true indeed since, for any m G Z 



g -(^+2,+2)/4 g (n-l)(n-2)/4 (g^ g)°o / PY^ ^ , ^-1 . ?> ^-2^2 



m /x m— 1 

y^(_!)i g -i(m-j)/2 xj = ( -(m-l)/2 x ^ = TT ft _ (m -l)/2 + ^ 

£Z (g;</)m-i(g;ff)i A = A 

This concludes the proof. □ 

5.2 A bilateral second order difference equation 

We know that the sequence u n = ) u+n (2w;q) obeys (|46)) . Applying the substitution 
q~ v ~ x = z, w = q% + if3, one finds that the sequence 

v n = q- n/A u n = q- n/A )„ + n{2q {2 » +1)/A P;q) (55) 

(?;?) c 

fulfills 

q<. n -Wpv n + (q n - z) v n+l + q n/2 f3v n+2 = 0, n G Z. (56) 
Remark 23. One can as well consider the unilateral second order difference equation 

(1 - z)v x + f3v 2 = 0, g ( "" 1)/2 /3t; n + (q n - z) v n+1 + q n/2 (3v n+2 = 0, n = 1, 2, 3, . . . . 

From (152]) it can be seen that the sequence {v n } given in ([35]) is square summable over 
N. Considering the Wronskian one also concludes that any other linearly independent 
solution of ( 15 6 p cannot be bounded on any neighborhood of +oo. Hence the sequence 
v n , iieN, solves the eigenvalue problem in £ 2 {N) iff vq = 0, i.e. iff ) v {2w\ q) = 0. In 
terms of the new parameters 0, z this condition becomes the characteristic equation 
for an eigenvalue z, 

(z' 1 ; q)oo <M; z -1 ; q, -z~ 2 /3 2 ) = 0. 
This example has already been treated in [TBI Sec. 4.1]. 
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For the bilateral equation it may be more convenient shifting the index by 1 in 
fl56|) . This is to say that we are going to solve the equation 

q^'^v^ + (q n -z)v n + q n/2 l3v n+l = 0, n G Z, (57) 

rather than (156|) . Denote again by J = J((3,q), with /3 G M and < q < 1, the 
corresponding matrix operator in £ 2 (Z). One knows, however, that J(—j3,q) and 
J(/3, g) are unitarily equivalent and so, if convenient, one can consider just the values 
f3 > 0. In equation (}5Tj) . z is playing the role of a spectral parameter. Using a notation 
analogous to flB]) (now for the bilateral case), this means that 

w n = q n/2 (3, X n = q n , and ( n ■= z - X n = z - q n , n G Z. (58) 

Notice that for a sequence {7 n } obeying 7„7 n+ i = u> ra , Vn G Z, one can take 

2 _ fc-1 2 _ kn2 

72/c-i — ? , 72fc — q P ■ 

Since the sequence {w n /(X n + 1)} is summable over Z, the Weyl theorem tells us 
that the essential spectrum of the self-adjoint operator J(/3, q) contains just one point, 
namely 0. Hence all nonzero spectral points are eigenvalues. 

Proposition 24. For < q < 1 and j3 > 0, the spectrum of the Jacobi matrix 
operator J((3,q) in £ 2 (1*), as introduced above (see §EE) ), is pure point, all eigenvalues 
are simple and 

spec p J(/3,g)= (-/3 2 g z +)Ug z 
Eigenvectors Vm^ corresponding to the eigenvalues q m , m G 7L, can be chosen asVm^ = 

= 1 (m - k)/4 )- m+k (2q- {2m+1)/4 P; q)- 
They are normalized as follows: 



IK +) || 2 = £ ^ /2 ^(2g- (2m+1)/4 /3;g) 2 = (-^ m /3 2 ;g)oo, Vm g z. 

fc=— oo 

Eigenvector Vm corresponding to the eigenvalues —f3 2 q m , m G Z +; can fre chosen as 
.(-) _ /„,(-) 



/_i \fc „fc(fc-4m-l)/4 

„(-) - 1 L) q R- k (-a- m+k+1 B~ 2 - a) 



x 



001 (; -q- m+k+1 p- 2 ; q, -g-^+^/T 2 ) . (59) 



Remark 25. An expression for the norms of vectors v m ^ can be found, too, 



(-)||2 _ ( i\m -m(3m+l)/2 (~QP 2 '■, l) oo {^Q ™ ft ^ g) oo (~ g m+1 /3 2 j g) oo 



-g/3 2 ;g)oo(g m+1 ;g) 



o 0i(;-^ 2 ;g,-g-^/3 2 ) g Z h 
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But the formula is rather cumbersome and its derivation somewhat lengthy and this 
is why we did not include it in the proposition and omit its proof. 

Proof. We use the substitution z = q~ v where v is in general complex. The right hand 
sides in (jl2j) can be evaluated using and (H2|) . Applying some easy simplifications 
one gets two solutions of fl57|) : 

v n = q-^ n ^) v+n {2q^l^-q), v n = (-l)V ( ^ +n)/4 )- v - n (2q^l^- q), n G Z. 

One can argue that in the bilateral case, too, all eigenvalues of J(f3,q) are simple. 
In fact, the solution {v n } asymptotically behaves as 

(9;9)oo 

For any other independent solution {y n } of 05 7p . q n l 2 (y n v n +i — Vn+iVn) is a nonzero 
constant. Obviously, such a sequence {y n } cannot be bounded on any neighborhood 
of +00. A similar argument applies to the solution {v n } for n large but negative. In 
particular, one concludes that z = q~ v is an eigenvalue of J(/3, q) if and only if {v n } 
and {v n } are linearly dependent. 

Using ([50]) one can derive a formula for the Wronskian, 

W (v, v) = q k,2 (3 (v k v k+1 - v k+ iv k ) 

= (-l) fc+1 /3g-^ 1 )/ 4 (j, +fe (2g( 2 - 1 )/ 4 /3; q)^ M (2q^^; q) 

+ j I , +fc+1 (2g^- 1 )/ 4 /3;g)j_^(2g^- 1 )/ 4 /3;g)) 

Thus z is an eigenvalue if and only if either q) 00 (qz; q) ^ = or (— z~ l (3 2 ; q) x = 
0. In the former case z G q z , in the latter case — z G f3 2 q z+ . 

Thus in the case of positive eigenvalues one can put v = —m, with m G Z. With 
this choice, {v k } coincides with {f^}. Notice that then the linear dependence of the 
sequences {v k } and {v k } is also obvious from (H3|) . Normalization of the eigenvectors 
Vm^ is a consequence of f )54p . 

As far as the negative spectrum is concerned, one can put, for example, r = 
— (irr + log/3 2 )/ logg and v = r — m, m G Z + . Then the sequence 

v k = q-^- m+k ^) T - m+k {-2iq^ 2m+ V ]q ), k G Z, 

represents an eigenvector corresponding to the eigenvalue — f3 2 q m . But it is readily 
seen to be proportional to the RHS of (15T?|) whose advantage is to be manifestly real. 

Finally let is show that can never be an eigenvalue of J((3,q). We still assume 
> 0. For z = 0, one can find two mutually complex conjugate solutions of (1571) 
explicitly. Let us call them v± iU , n G Z, where 
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Clearly, 

v±, n = i ±n q~ n/4 (1 + (q n/2 )) as n -+ +00. 
Using the asymptotic expansion one can evaluate the Wronskian getting 

W (v+, v.) = q n/2 (3 (u +>n u_, n+ i - v + , n+1 v_, n ) = -2iq- l ' A l3. 

Hence the two solutions are linearly independent. It is also obvious from the asymp- 
totic expansion that no nontrivial linear combination of these solutions can be square 
summable. Hence cannot be an eigenvalue of J(j3, q) whatever j3 is, and this con- 
cludes the proof. □ 

So one observes that the positive part of the spectrum of J(/3, q) is stable and does 
not depend on the parameter (3. This behavior is very similar to what one knows 
from the non-deformed case. On the other hand, there is an essentially new feature 
in the q-case when a negative part of the spectrum emerges for j3 7^ 0, and it is even 
infinite-dimensional though it shrinks to zero with the rate j3 2 as f3 tends to 0. 



6 Q-confluent hypergeometric functions 

In this section we deal with the q-confluent hypergeometric function 

( a ;<7)fc k 

^{a-O-q.q-z) = > l-irg-" „ \ , ^z ■ 



k=0 

It can readily be checked to obey the recurrence rules 

a+7/1 _ 7-q+1\ / 7 

y V 1 y ' ., J. r „a. „<Y+2. „ „<y+2 \ / 1 H z 



(1 - g7) (1 - g7+i) 

+ i0i(g a ;g 7 ;g,g^) = O (60) 
and 

i<Mg^ +1 ; g 2 " 7 ; g, *) + 9 (? ~ + ^ nM?-*- 1 ; <T 7 ; ?, z) 

q(q- q 1 - g 1 " 7 + 1) g - g 7 \ x 

+ « z 1^1(9 ;y 7 ;g,^) = 0. 



Put, for n G Z, 

^ = (g" +7 ; g)oo i<M<f ; g" +7 ; 9, -^ +7 ^), (61) 

(„n+7-a. ~\ 

^ = g-a(n+7)-(n+7-l)(n+7-2)/2 ^^a-n-7+1. ? 2-n- 7; g? _ qz y ^ 

Here z, a, 7 G C, g 7 ^ g z . The recurrence rules imply that both {</?„} and {tp n } solve 
the three-term difference equation 

q a +1+ n-l ^ _ q -y-a+nj _ ^ _ ?7+ „-l + ? 7+n-l^ ^ + ^ — Q, n G Z. (63) 
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Lemma 26. The sequences {tp n } and {VVi} defined in l[61\) and HM) . respectively, 
fulfill 

Poipi ' Vi4>o = q- a ^ +1) -^ (j - 1 \q J - a+1 ;q)oo(-q a z;q) QO z-\ (64) 
Alternatively, [64\ ) can be rewritten as 

Mq a ; g 7 ; q, g 7- ^) i<Mg a ~ 7 ; g 1-7 ; q, g 1- ^) 

+ (i-ql^) (i - lMqa; g7+1; g ' qJ ~ a+lz) ^i(<i a - 7+1 ; <? 2 " 7 ; ?. q 1 '^) = (*; 9)00. 

Proof. Checking the Wronskian of the solutions and ^>„ one finds that 

? i n(n _ 1)+(Q!+7)n (?7 _ Q+1 . ? ^ _ = c (65) 

is a constant independent of n. In particular, v^oV'i — { Pi' 1 Pq — C. It is straightforward 
to examine the asymptotic behavior for large n of the solutions in question getting 
(p n = 1 + 0(q n ) and 

^ n = g -|n(n-l)-( a + 7 -l)n-i( 7 -l)( 7 -2)-a 7 ^l- 7 -n ^ ^ + _ 

Sending n to infinity in fl65l) one finds that C equals the RHS of (I64j) . □ 
Proposition 27. For a, 7, 2; G C, 

? — 7^ — zv^eztW = 771 — 1x3773 i<Mg ; g 7 ; g, -gV)- 



( g7 -a+fe+i. g2)oo (1 - (1 - 2)qry+*-i) / ((1 - z) g 7 ; g ) c 

fc=i/ 



(66) 



Proof. The both sides of the identity are regarded as meromorphic functions in z. 
Setting I1117 to a constant, the both sides tend to 1 as Re 7 tends to +00. In virtue 
of Lemma El it suffices to verify that the sequence 

(fff+n-l. n \ 



1 (j )l (q a ;q^ n - l ;q,-q^ n - 1 z), n G N, 



((l-z)g7+n-i ;g ) c 

satisfies the three-term recurrence relation F n — F n+ i + s n zF n+2 = 0, n G N, where 

Sn ~ (1 - (1 - z)q<+ n - r ) (1 - (1 - z)q< +n ) ' 

Since 7 here is arbitrary one can consider just the equality for n — 1. But then the 
three-term recurrence coincides with (1601) (provided z is replaced by —z). □ 

Let us now focus on equation f )63|) . One can extract from it a solvable eigenvalue 
problem for a Jacobi matrix obeying the convergence condition (I14p . 
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Proposition 28. For a G R and 7 > — 1, Ze£ J = J(cr, 7) 6e i/ie Jacobi matrix operator 
in £ 2 (N) defined by and 

n>„ = ~ sinh^^-^Vl-?^ 7 . ^ = g™" 1 . (67) 
Tnen 2 7^ is an eigenvalue of J (a, 7) i/ and on/y i/ 

(cosh 2 {a/2)z~ 1 50)00 i0i(g~' r cosh 2 ((7/2)z~ 1 ;cosh 2 ((T/2)z" 1 ;a,-sinh 2 ((7/2)^ 1 ) = 0. 
Moreover, if z 7^ solves this characteristic equation then the sequence {v n }^ =l , with 

V n = q -hn + Hn- 3) Smh>)(2,)- f cQsh2 ^ N 

x ^(q^ cosh 2 (|) z- 1 ;?" cosh 2 (|) z" 1 ; g, -g" sinh 2 (£) z' 1 ) , (68) 

a corresponding eigenvector. 
Remark 29. Notice that the matrix operator J{a, 7) is compact (even trace class). 
Proof. First, apply in (1631 the substitution 

7 = 7 + a, z = q 13 , u n = q~ an u n , 
and then forget about the tilde over 7 and u. Next use the substitution 

q p/2 = tanh(0, q a = g~ 7 cosh 2 Q z~ x , u n = (j) n u n . 
where {<f) n } is a sequence obeying 

^ n = g (/3+7+n-l)/2 / 1 _ g7 +„ _ 
0n+l 

Up to a constant multiplier, n 2 = q-P n -"/ n -2 n ( n ~ 3 ) (qY+ n ^ q)^. We again forget about 
the tildes over z and u, and restrict the values of the index n to natural numbers. If 
u = then the transformed sequence {v,k}T = i solves the Jacobi eigenvalue problem 
(JHJ) with w n and A n given in (l6Tj) . 

Further apply the same sequence of transformations to the solution (p n in (l6Tj) . 
Let us call the resulting sequence {v n }. A straightforward computation yields (|68|) . 
Clearly, the sequence {t^; k > 1} is square summable. On general grounds, since 
J(a, 7) falls into the limit point case, any other linearly independent solution of the 
recurrence in question, (ITT)) , cannot be square summable. Hence the characteristic 
equation for this eigenvalue problem reads no = 0. This shows the proposition. □ 

Remark 30. In the particular case 7 = the characteristic equation simplifies to the 
form 

(coshV/2)*- 1 ; q) oo (- sinh V/2K 1 ; q) oo = 0. 
Hence in that case, apart of z — 0, one knows the point spectrum fully explicitly, 

specJ(a,0)\{0} = {q k cosh 2 (a/2); k = 0, 1, 2, . . .} U {-q k sinhV/2); k = 0,1,2,...}. 



26 



Remark 31. Of course, Proposition 1281 can be as well derived using formulas ( JT6 
(I18j) . while knowing that (114")) is fulfilled. To evaluate C, n (z) one can make use of ()66 
Applying the same series of substitutions as above to equation f l66|) one gets 



g2( fc -7)-| S inh(a) (g 7+fc ;g 2 ), 
2(g7+fc+ i. g2)oo(gfc -i_ ;z) 

(cosh 2 (a/2) z-\q) 



1 OO 



fc=l 

-1 



(^ -1 ;<7)oo 

Then a straightforward computation yields 



i0i (q 7 cosh 2 (t^jz cosh 2 (0 ^ g, - sinh2 (^) 21 



Cn(«) = — • w x 7 — — ^ = o, 1, 2 

sinh(o-) (z 1 ;q) 00 



with f n being given in 



Acknowledgments 

One of the authors (F.S.) wishes to acknowledge gratefully partial support from grant 
No. SGS12/198/OHK4/3T/14 of the Grant Agency of the Czech Technical University 
in Prague. 



References 

M. Horowitz, I. A. Stegun: Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, (Dover Publications, New York, 1972). 

N. I. Akhiezer: The Classical Moment Problem and Some Related Questions in 
Analysis, (Oliver & Boyd, Edinburgh, 1965). 

B. Beckerman: Complex Jacobi matrices, J. Comput. Appl. Math. 127 (2001) 
17-65. 

J. Gard, E. Zakrajsek: Method for evaluation of zeros of Bessel functions, J. Inst. 
Math. Appl. 11 (1973) 57-72. 

G. Gasper, M. Rahman: Basic Hypergeometric Series (Cambridge University 
Press, Cambridge, 1990). 

W. Gautschi: Computational aspects of three-term recurrence relations, SIAM 
Rev. 9 (1967) 24-82. 

I. S. Gradshteyn, I. M. Ryzhik: Table of Integrals, Series, and Products, Edited 
by A. Jeffrey and D. Zwillinger (Academic Press, Amsterdam, 2007). 



27 



Y. Ikebe: The zeros of regular Coulomb wave functions, Math. Comp. 29 (1975) 
878-887. 

Y. Ikebe, Y. Kikuchi, I. Fujishiro: Computing zeros and orders of Bessel functions, 
J. Comput. Appl. Math. 38 (1991) 169-184. 

Y. Ikebe, Y. Kikuchi, I. Fujishiro, N. Asai, K. Takanashi, M. Harada: The eigen- 
value problem for infinite compact complex symmetric matrices with application 
to the numerical computation of complex zeros of Jq(z) — iJ\{z) and of Bessel 
functions J m (z) of any real order m, Linear Alg. Appl. 194 (1993) 35-70. 

F. H. Jackson: On generalized functions of Legendre and Bessel, Trans. Roy. Soc. 
Edinburgh 41 (1903) 1-28. 

Y. Miyazaki, Y. Kikuchi, DongSheng Cai, Y. Ikebe: Error analysis for the com- 
putation of zeros of regular Coulomb wave function and its first derivative, Math. 
Comp. 70 (2001) 1195-1204. 

M. Rahman: An addition theorem and some product formulas for q-Bessel func- 
tions, Can. J. Math. 40 (1988) 1203-1221. 

F. Stampach, P. St'ovfcek: On the eigenvalue problem for a particular class of 
finite Jacobi matrices, Linear Alg. Appl. 434 (2011) 1336-1353. 

F. Stampach, P. St'ovicek: The characteristic function for Jacobi matrices with 
applications, arXiv: 1201 . 1743^ 1 [math.SP] (submitted). 



28 



